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ABSTRACT 

We report the results of Australia Telescope Compact Array (ATCA) obser- 
vations of the Westerlund 1 (Wdl) region in the SiO v — 1, J — 1-0 and H 2 
6i6~523 maser lines, and we also report the analysis of maser properties of red 
supergiants (RSGs) associated with 6 massive clusters including Wdl. The pri- 
mary purpose of this research is to explore possibilities of using maser emission 
for investigating the nature of massive clusters and associated RSGs. The SiO 
v = 1, J = 1-0 and H2O 6i6-523 maser lines are detected toward 2 of 4 known 
RSGs in Wdl. The large velocity ranges of maser emission are consistent with 
the RSG status. RSGs with maser emission tend to exhibit redder log(F 2 i/F 12 ) 
and [K-12.13] colors compared to RSGs with no maser emission. The mass-loss 
rates derived from dust radiative transfer modeling suggest that RSGs with maser 
emission tend to exhibit larger mass-loss rates compared to RSGs with no maser 
emission. In an extended sample of 57 RSGs in 6 massive clusters, detections in 
the SiO line tend to homogeneously distribute in absolute luminosity L, whereas 
those in the H2O line tend to distribute in a region with large L values. 



Subject headings: masers — (Galaxy:) open clusters and associations: individual 
(Westerlund 1) — stars: late-type — stars: mass- loss — (stars:) supergiants - 
radio lines: stars 
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Introduction 



Until now, about 200 red supergiants (RSGs) have been identified in the Galaxy (ILevesque 



20101 ) . but a majority of known RSGs are isolated field stars, which pose difficulties in ob- 



taining critical stellar parameters, such as distance, age and luminosity. However, in the 



Westerlund 


1 (Wdl; Menee 


& Tacconi- 


Garman 2007) 


, Red Supergiant Cluster 1 (RSGC1; 


Fieer et al. 


2ooe 


), RSGC2 ( 


Davies et a 


2007 


), RSGC3 ( 


Alexander et a 


2009: Clark et al. 


2009 


), RSGC4 ( 


Neeueruela et al. 


2010a) 


, and RSGC5 


(Negueruela et al. 


2011). These clus- 



ters that harbor multiple RSGs at a time are suitable sites for investigating the nature of 
RSGs, because RSGs in a single cluster share a coeval condition and the same distance. 

It is notable that RSGs often exhibit maser emission of the SiO, H2O and OH molecules. 
Since at the frequencies of these maser lines, interstellar extinction is negligible, maser emis- 
sion potentially could be a useful tool for investigating RSGs hidden in the galactic plane. 
However, maser emission of RSGs embedded in massive clusters still has not been well- 
studied. In order to explore possibilities of using maser emission for investigating massive 
clusters and embedded RSGs, we have made maser searches toward 5 massive clusters, 
which are observable from the northern hemisphere, in the SiO (v = 1 and 2 , J = 1-0) and 



H 2 (22 GHz) maser lines (INakashima fc Deguchill2006l ; iDeguchi et al.ll2010l ). and until now 
roughly a dozen of maser sources of RSGs have been found in the 5 massive stellar clusters. 
In addition, the results of an Si O maser survey of Ga lactic RSGs in the SiO v — 1, J = 2-1 
line has been recently reported (jVerheyen et al.l 120121 ). 



In this paper, we report the results of interferometric maser observations using the 
Australia Telescope Compact Array (ATCA) toward the massive cluster Wdl in the SiO 
and H2O maser lines at 43 GHz and 22 GHz, respectively. The purposes of this research 
are (1) to search SiO and H2O maser emission toward the Wdl region and (2) to explore 
possibilities of using maser emission for investigating the nature of massive clusters and 
associated RSGs. To clarify the maser properties of RSGs in massive clusters, we analyze 
the previous data of maser observations of massive clusters, in conjunction with the present 
observation of Wdl. 

The structure of this paper is as follows. In Section 2, we describe the details of the 
observations of Wdl. In Section 3, we summarize the results of SiO and H 2 maser obser- 
vations of Wdl. In Section 4, we discuss the velocity information derived by the present 
observations. In section 5, we discuss the infrared properties of RSGs in Wdl and other 
massive clusters. Finally, the main results of the present research are summarized in Section 
6. 
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2. Observations and Data Reductions 

The radio interfero metric observation in the SiO v — 1, J — 1-0 line towards the 
Wdl region was made with ATCA on 2008 November 20 and 21 (pr oject code: C1796, PI : 



Nakashima, J.). Since the angular size of the Wdl region (2' x 4'; iBorgman et al.l Il970l ) 
is much larger than that of a single field-of-view (FOV) of an ATCA 22 m antenna at the 
observing frequency of 43 GHz (HPBW~ 66"), we observed 8 FOVs with different center 
positions to cover all RSGs and RSG c andidates in Wdl (center po s itions of observed FOVs 



are summarized in Table 1). In Wdl, iMengel fc Tacconi-Garmanl (120071 ) have identified 4 



RSGs (W 20, W 75, W 26 and W 237) through their near infrared spectroscopy, and all 
the 4 known RSGs were covered in the present observation (coordinate values of the known 
RSGs are summarized in Table 2). In addition to 4 known RSGs, to make doubly sure, a 
majority of i^-band bright stars in the Wdl region (i.e., possible mass- losing evolved stars) 
are covered with 8 FOVs. In Figure 1, we present the locations of 8 observed FOVs (see, the 
blue circles in Figure 1) and 4 known RSGs. 

The weather condition was stable throughout the entire observation run. The observa- 
tion used the EW352/367 array configuration consisting of six 22 m antennas. The baseline 
lengths ranged from 50 to 4400 m. The rest frequency of the SiO v = 1, J = 1-0 line 
(43.122090 GHz) was centered at a correlator window with a bandwidth of 32 MHz. The 
frequency coverage was ranged from 43.148 to 43.116 GHz. The 256 channel correlator gave 
the effective velocity coverage and resolution of 222.6 km s _1 and 0.869 km s _1 , respectively. 
The quasar 1253—055 (3C279) was observed at the beginning of the observation for band- 
pass and flux calibrations. A nearby strong SiO maser source, IRAS 16105—4205 was also 
observed at the beginning for checking reliability of the system (the SiO v — 1, J — 1-0 line 
is detected for the first ti me toward this obje ct, even though the strong v — 1, J — 2-1 line 



was previously detected; lHaikala et al.lll994 ; see, Appendix A, for details). The resulting 



uncertainty of the absolute flux density is roughly within 15-20%. Observations were made 
in a snap-shot mode, which observed 8 FOVs in turn. The observation was periodically 
interleaved with the nearby gain calibrators, 1646—50, to track the phase variations over 
time (the phase calibrator was observed for 3 minutes every roughly 30 minutes). After all 
positions were observed, the array slewed to the first position to loop again. At least 3 loops 
were repeated for each position throughout the whole observing session. The on-source in- 
tegration time at each FOV was 10 minutes in each single loop except for positions 2 and 3, 
which took 15 minutes each (this is because positions 2 and 3 are the most crowded region in 
Wdl, a better signal-to- noise ratio helps to avoid confusion in source identification). Total 
integration time on position 1, positions 2-3, and positions 4-8 were 0.7, 1.0, and 0.5 hours, 
respectively. And anot her 1.6 hours we re used for calibrations. Visibility data were cali- 



brated using MIRIAD ( iSault et al.lll995l ). following a standard calibration procedure. Then 
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we applied a self-calibration technique using the MIRIAD task, self cat to enhance the signal- 
to-noise ratio. Image processing of the data was also performed with MIRIAD by applying 
robust weighting, which compromises between the natural and uniform weighting (we ap- 
plied "robust=0.5"). As we observed different FOVs in a snap-shot mode, the u-v coverages 
are different from FOV to FOV. Therefore, synthesized beams at each FOV have slightly dif- 
ferent sizes and position angles; the information about synthesized beam sizes is summarized 
in Table 1 (note: averaged synthesized beam size is 6.5" x 2.4"). We checked the detection of 
the continuum emission by integrating over all emission-free channels (the integration range 
is about 25 MHz), but did not detect it. The upper limit is 1.5 x 1CT 3 Jy beam -1 . 

Archival ATCA data of a previous observation in the H2O 6i6~523 line toward Wdl 
(project code: C1619, PI : Sean Doughertjj^) were used for the present analysis. The H 2 
maser observation was made from 2006 December 6 to December 10 using the 6A array 
configuration. The six antennas form baselines ranged from 340 to 5900 m. The observed 
H2O line at 22.2351204 GHz is centered at a correlator window with a bandwidth of 16 MHz. 
The frequency coverage is ranged from 22.243 GHz to 22.227 GHz. The 512 channel corre- 
lator gave an effective velocity coverage and resolution of 215.7 km s _1 and 0.421 km s _1 , 
respectively. The HPBW of a single antenna is 126" at the observing frequency. With 
this beam size of a single antenna, almost the entire cluster region was covered only with 2 
FOVs. We present the locations of observed FOVs in Figure 1 (see, the green circles) and 
Table 1. The observation was made in a snap-shot mode as well as the SiO maser observa- 
tion. A nearby phase calibrator, 1646—50 was observed every 10 minutes. A flux calibrator 
used was 1253—055 (3C279). Calibration and imaging of the data were performed using 
MIRIAD following the same methodology as adopted for the SiO data. We applied robust 
weighting (robust=0.5) giving an averaged synthesized beam size of 1.3" x 0.4". We did not 
detect any continuum emission (integration range is roughly 10 MHz); the upper limit is 
2.4 x 10~ 3 Jy beam -1 . 

3. Results of Maser Observations of Wester lund 1 

3.1. Sources Identification 

First of all, we briefly describe how we identified stellar maser sources in the data cubes 
of the present observation. Source identification processes are common in both SiO and 



1 The OH 1612 MHz line was also observed in this project, and the emission seems to be detected toward 
Wdl according to our brief inspection. However, as the data quality is not good enough for further analysis, 
we do not use it in the present research. 
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H2O observations. We first sorted out all emission features detected above a 3 a level in 
each velocity channel of calibrated data cubes, and obtained the coordinate values of each 
emission feature by fitting a two-dimensional Gaussian function using a MIRIAD task imfit. 
At this stage, roughly 10 emission sources are selected as possible detections in each data 
set of the SiO and H2O observations. Then, we excluded artificial features in the following 
procedure. 

Generally speaking, true maser lines exhibit a linewidth larger than, at least, 1-2 km s" 1 , 
which is equivalent to 2 and 3 channels at 43 GHz and 22 GHz, respectively (velocity resolu- 
tions of the present observation at 43 GHz and 22 GHz are 0.869 km s _1 and 0.421 km s" 1 , 
respectively, as we mentioned in Section 2). Therefore, we looked for the emission features 
that consistently appear at the same position within uncertainty of a synthesized beam size 
and survive respectively in more than 2 and 3 consecutive velocity channels at 43 GHz and 
22 GHz, so that we could exclude spiky features with a very narrow linewidth, which are 
most likely artificial. 

Then, we also doubly checked the reliability of detections by cross-checking with 2MASS 
images. Since the photosphere of RSGs exhibits a low effective temperature of about 3000 K, 
a bright near IR counterpart with a red IR color should be found at the positions of maser 
sources if emission features are really originated in a RSG. We cross-checked between our 
radio images and 2MASS /f-band images to fi nd IR counterparts. A s the angular resolution 



of 2MASS images is about 2.0" in if-band (ISkrutskie et al.l 120061 ). the resolution of the 



present observation seems to be sufficient enough for positional comparison. In considerations 
of beam sizes of both 2MASS and the present observations, we defined that emission is real 
if the 2MASS point sources with a red color are lying within 4.5" and 2.0" with respect to 
the radio positions of emission features at 43 GHz and 22 GHz, respectively. Through these 
steps, finally we identified two real emission sources detected each in the SiO and H 2 lines 
(the SiO and H 2 lines are detected toward the same two sources). The details of the results 
are summarized in Section 3.2. 



3.2. Results of the SiO and H 2 Maser Observations 

As mentioned in Section 3.1, we finally detected the SiO v = 1, J = 1-0 line toward 2 of 
4 known RSGs (W 26 and W 237) in Wdl, and also detected the H 2 6i 6 -5 23 line toward the 
same two RSGs. In addition, we inspected across entire observing regions to doubly check 
further maser detections of unknown mass-losing evolved stars, but no other detections were 
found at a 3 o level. In Table 3 we present the coordinate values of detected maser sources, 
peak and velocity integrated intensities, intensity peak velocities and rms noise levels. In 
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Table 4 we present rms noise levels toward two known RSGs with no maser emission. In 
Figure 2 we present the line profiles of detected maser lines. 

Both W 26 and W 237 exhibit asymmetric profiles (asymmetric with respect to the 
peak velocity) with a wide velocity range; the velocity ranges of W 26 and W 237 are 
roughly 40 km s _1 and 23 km s _1 , respectively. Generally speaking, RSGs exhibit a wider 
velocity range than that of asymptotic giant branch (AGB) st ars, which exhibit a velocity 
range of, at most, 10 km s _1 in the SiO maser lines (see, e.g., iNakashima fc Deguchil 120071 ; 
Zhang et al.ll2012al ). Therefore, the velocity range of SiO maser emission is consistent with 
the RSG status. Additionally, these two stars clearly exhibit asymmetric profiles (note: the 
black arrows in Figure 2 indicate weak co mponent above a 5 a level). This is also typica l 



characteristics often seen in RSGs (see, e.g.. INakashima &: Deguchill2007l ; IZhang et al.ll2012al ). 



Both W 26 and W 237 exhibit clearly different intensities in the SiO line (the peak intensities 
of W 26 and W 237 are 0.12 Jy and 1.54 Jy, respectively), while both stars exhibit almost 
the same intensity in the H2O line (the peak intensities of W 26 and W 237 are 0.73 Jy and 
0.72 Jy, respectively). Since the two stars are at the same distance as in the same cluster, 
the difference in intensity directly reflects the difference of the absolute intensity of maser 
emission. This difference in intensity in the Si O maser line may be explained by the tim e 
variation in intensity of maser emission (see, e.g.. lKamohara et al.ll2005l ; IZhang et al.ll2012al ). 



We shall note that the intensity of the H2O maser line is somewhat weaker than that 
of three well-known nearby (isolated) RSGs: i.e., VY CMa (1.2 kpc), VX Sgr (1.6 kpc) and 
NML Cyg (1.6 kpc) exhibit the pe ak intensities of 199.8 Jy, 167.2 Jy and 18.8 Jy (these 



flux v alues are scaled at 3.5 kpc, iPashchenko et al.l 120061 ; IChoi et al.l l2008t IZhang et al. 



2012bh . while W 237 and W 26 exhibit 0.72 Jy and 0.73 Jy. Here, we find roughly two 
order of differences in intensity between the two cases. Similarly, the SiO maser emission 
(v = 1, J = 1-0) of the three well-known nearby RSGs exhibit systematically larger fluxes 
than RSGs in Wdl: VY CMa, VX S gr and NML Cyg exhibit t he peak intensities of 5.4 Jy. 



31.4 Jy and 2.6 Jy (scaled at 3.5 kpc. lBoboltz &: Claussenll2004t IZhang et al.ll2012al ; ISu et al. 



20121 ). while W 237 and W 26 exhibit 1.54 Jy and 0.12 Jy, respectively. However, these 
differences in intensity between nearby RSGs and RSGs associated with Wdl are partly 
due to the uncertainty of distances. In fact, the three representative nearby RSGs exhibit 
relatively large 12 /xm fluxes as well as the large maser intensities: i.e., the 12 /im fluxes of 
VY CMa, VX Sgr and NML Cyg are 1,166 Jy (IRAS), 433 Jy (MSX) and 587 Jy (MSX; 
these flux values are scaled at 3.5 kpc), while those of W 237 and W 26 are 101 Jy, and 
326 Jy, respectively. The systematically large mid-infrared fluxes of the nearby RSGs means 
that the distances used for scaling may include a large uncertainty (even though a large 
dispersion in infrared intensities should be also taken into account). Although the distances 
to the nearby RSGs are based on the trigonometric parallax measurements giving a reliable 



-7- 



distances, as we discuss later in Section 4, the distance to Wdl possibly includes a large 
uncertainty. Therefore, for more precise discussions on absolute maser intensities, we need 
more reliable distance information about Wdl. 

On the other hand, the intensity ratios of the H2O to SiO maser intensities are distance- 
free parameters, and the nearby RSGs seem to exhibit a relatively large value on average: the 
peak intensity ratios of the H 2 to SiO maser lines of VY CMa, VX Sgr and NML Cyg are 
36.9, 5.3 and 7.3, while those of W 237 and W 26 are 0.47 and 6.1. Even though the results 
is not very obvious (because only VY CMa and W 237 exhibit extreme values), one may 
consider that this relatively weak intensity of the SiO maser line of RSGs in Wdl is somehow 
originated in the cluster environment; for example, the outer layer of the envelope may be 
disturbed by interactions with ambient materials in the condensed cluster environment. 

It is notable that SiO maser emission of W 26 and W 237 exhibits larger velocity ranges 
than H2O emission. In the case of spherical envelopes, the velocity range of SiO emission 
theoretically should be smaller than that of H2O emission, because molecular gas is usually 
still being accelerated in the radial direction in the innermost region (2-3 R*), from which SiO 
maser is emitted. Since maser emission cannot be amplified in the direction of acceleration, 
SiO masers are usually amplified only in the tangential direction to the radial direction. 
Therefore, in the case of a spherical envelope, intensity peaks of SiO maser emission are 
found at velocities close to the systemic velocity. On the contrary, at the H2O maser region 
(20-30 R*) the motion of molecular gases usually already reaches to almost the terminal 
velocity and, in addition, gas is usually turbulent. Therefore, the maser can be amplified 
in any direction. As a result, the velocity range of H2O maser emission is generally wider 
than that of the SiO maser emission. The present result possibly suggests that circumstellar 
envelopes of observed RSGs somewhat deviates from spherical symmetry. We shall note, 
however, that this discussion may be affected by the time variation of the line profile and 
sensitivity of the observation if weak components exist below the detection limit (for example, 
as mentioned later in Section 4, the convective motion could alter the line profile). 



4. Systemic Velocities of Maser Sources 

In the case of AGB stars, it is well known that the peak velocity of SiO m aser emission 



repre sents the stellar systemic velocity within an uncertainty of 2-3 km s 1 ( 1 Jewell et al. 



1984 ). In the case of RSGs, however, the peak intensity clearly does not represent the 



systemic velocity of the stars, because the line profiles are not symmetric in velocity. For 
that reason, here we estimate the systemic velocities of RSGs in Wdl by assuming that 
the median velocity of a line profile is close to the stellar systemic velocity: specifically, we 
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calculated the average of the highest and lowest velocities of 5 a channels. In the case of 
SiO maser emission, the derived values of the median velocities for W 26 and W 237 are 
—38.1 km s _1 and —41.7 km s _1 , respectively. Similarly, the median velocities obtained 
from the H 2 maser profiles are —39.3 km s _1 and —48.5 km s _1 for W 26 and W 237, 
respectively (since W 26 exhibits a single-peak profile in the H2O line, here we simply give the 
peak velocity). These derived values are not exactly consistent with the systemic velocities 
obtained by infrared high- dispersion spectroscopy of the CO band-head at 2.29 jttm (near- 
infrared velocities are —44.5 km s _1 and —57.5 km s _1 for W 26 and W 237, respectively; 



Mengel fc Tacconi-Garmanll2009[ ). The differences between the present observation and near- 
infrared measurement are roughly 5-6 km s _1 and 9-16 km s _1 for W 26 and W 237, 
respectively. 

As we mentioned in Section 2, we checked the reliability of our velocity calibration by ob- 
serving a known maser source with a known systemic velocity, and c onfirmed that the pr esent 



velocity measurement is consistent with a previous measureme nt (IHaikala et al. 



1994). On 



the o ther hand, the resolution of near-infrared observations (IMengel fc Tacconi-Garman 
20091 ) was also pretty high (R ~ 9000; this corresponds to Av = 33 km s _1 ), and the 
accuracy of the infrared velocity measurement must reach to, at least, a several km s _1 
level, because a series of absorption lines were simultaneously observed and the derived ve- 
locities were averaged out. Therefore, inconsistency in the velocities seems to be originated 
in astrophysical reasons. In fact, the atmosphere of RSGs are strongly affected by connec- 
tive motions: the maximum expansion ve locities of RSG envelo pes could be often beyond 
30-40 km s _1 due to convective motions ( IJosselin &: Plezl 120051 ). In addition, the infrared 
velocity determined by observing the CO band-head may be affected by the pulsation of 
RSGs; for example, in the case of M-ty pe miras, the pulsatio n is known to shift the position 
of the CO band- head up to 15 km s _1 ( iScholz fc Woodll2000l ). Presumably, the difference in 
velocities up to 16 km s _1 may be interpreted for these reasons. 

The average of two systemic velocities of SiO maser emission is calculated to be —39.9 km s" 
Similarly, in the case of H2O maser emission, the average is —43.9 km s _1 . These values are 
expected to be close to the systemic velocity of Wdl. Previously, two indirect measurements 
of the systemic velocity of Wdl have been made based on t he motion of ambient inters tellar 
gas (HI and CO) surr ounding Wdl: — 55^ km s -1 (HI; iKothes fc Dougherty! 120071 ) and 
-48 ± 14 km s" 1 (CO; iLuna et all 120091 ). The present values are not inconsistent with these 
values if we consider the uncertainty (the systemic velocities of Wdl are briefly summarized 
in Table 5). 

Using the systemic velocity derived, we calculated the kinematic distance to Wdl under 
the assumption of a flat circular rotation model of the Galaxy (here, we use the velocity 
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of SiO maser emission a s a representative value): in the calculation we followed the same 
methodology adopted by lKothes fe Dougherty! (120071 ). assuming the galactocentric distance 
of 7.6 kpc, and the derived kinematic distance is 3.0 ± 0.5 kpc (here we assume an uncer- 
tainty of 10 km s _1 in the velocit y of Wdl based on the above discussion). We also calculated 
the kinematic distance based on iReid et al.l (120091) for c ompa rison, and found no significant 
differences from those based on iKothes fc Dougherty ( 2007 ): the obtained kinematic dis- 
tances are summarized in Table 6. On the other hand, the obtained kinematic distance is 
slightly nearer than luminosity distances derived by independent methodolo gies using main 



sequence stars, OB supergiants and Woff-Rayet star s in Wdl (3.5-5.5 kpc, ICrowther et al. 



20061 : iBrandner et al.l 120081 ; iNegueruela et al.ll2010bl ). Therefore, in the following analysis, 



we adopt 3.5 kpc as the representative distance to Wdl, because this distance is not incon- 
sistent with both the kinematic and luminosity distances [at the distance of 3.5 kpc, absolute 
luminosities of W 26 and W 237 are respectively calculated to be 1.1 x 10 6 L Q and 2.3 x 10 5 
L based on the photometric data mentioned in Section 5.1; these values are roughly con- 
siste nt with a range of typical absolute luminosities o f RSGs (2 x 10 4 L — 6 x 10 6 L ; see, 
e.g.. Istothers k ChirJll999l : iMauron fc Josselmlboill )] . 



5. Discussion 

Our primary purpose of this project is to explore possibilities of using maser emis- 
sion for investigating massive clusters and associated RSGs. We previously discussed the 
kinematics of the cluster itse lf using the velocity information of RSGs with maser emission 



( iNakashima fc Deguchil 120061 ) . However, since the number of detection in the present obser- 
vation is not sufficient enough for discussing cluster kinematics, here we focus on the nature 
of RSGs. In the following subsections, firstly we compare infrared properties of 4 known 
RSGs to clarify the difference and similarity between detections and non-detections in the 
SiO and H2O maser observations. Secondary, we extend our analysis to RSGs in other mas- 
sive clusters, which previously have been observed in the SiO, H 2 and OH maser lines, so 
that we can more comprehensively consider the characteristics of RSGs with and without 



maser emission. 



5.1. Comparison between Detections and Non-Detections 



The infrared data used for the present analysis were ta ken from followin g data archives: 
Deep Near Infrared S urvey of the Southern Sky (DENIS; beul et al.l[l995l ). 2MASS Point 
Source Catalog (PSC; ISkrutskie et al.ll2006l ). Spitzer's Galactic Legacy Infrared Mid-Plane 



- 10 - 



Survey Extraordi naire (GLIMPSE; IWerner et al.ll2004l ). Midcourse Space Experiment (MSX; 
Egan et al.l 120031). Wide-field Infrared Survey Explorer (WI SE) All Sky Source Catalog 
rtWright et aJboid ) and AKARI/IRC Point Source Catalog Jlshihara et al.lboiof ). In con- 
junction with these archival data, we obtained 8-18 photometric data points per each RSG 
from 0.82 /im to 22 /im. The photometric data collected are given in Appendix B. 

It has been known that interstellar extinction toward Wdl is non-negligible in near- 
infrared wavelengths. Therefore, we adopted an extinction correction for p hotometric data 



below 10 /i m in the following way. F ir stly, we found Ej_ K values from iLevesque et al. 



( 120051 ) and iMengel fc Tacconi-Garmanl (120071 ) toward each RSG; their Ej_k values are 
based on intrinsic J-K colors determined by effective temperatures of theoretical models 
( ILevesque et al.ll2005l ; IMengel &: Tacconi-Garmanl 120071). Secondary, A K values toward each 
RSG are determined using the formula given in ICrowther et al.l (120061) [ note: some authors 
recently use a slightly different formula (see, e.g.. iNishiyama et al.l 120061 ). but the difference 
of the formulas does not affect the results of our discussion; see, some more details in Sec- 
tion 5.2]. Finally, der ived extinction at i f -ban d is extrapolated to other wavelengths using 
the formula given in llndebetouw et al.l ( 120051 ). The averaged extinction toward 4 known 
RSGs at if-band is 1.33. This averaged extinction is consistent with previous studies; for 
example, the extinction measur ements using main- seq uence and pre- main- sequence stars in 
Wdl shows Ak is roughly 1. 13 (jBrandner et al.l 120081 ). the result of Wolf-Rayet stars shows 



A K i s 1.01 (ICrowther et al.l 120081). and the results of OB supergiants are 1.20 (I Clark et al. 



20051 ) and 1.34 (INegueruela et al.ll2010bl ). Both original and interstellar extinction corrected 
spectral energy distribution (SED) diagrams are given in Figure 3. 

To estimate physical parameters of dust envelopes, SED profiles were fitted by one- 
dimensional dust radiative transfer models assuming centrally-heated spherical density dis- 
tributions using the DUSTY code ( llvezic et al.lll999l ). The central star (i.e., central heat 
source) was assumed to be a point source at the center of the SEDs and their spectral energy 
distributions were take n to be Planckian. We used the grain type of co ld (Sil-Oc) silicates 
dOssenkopf et al.lll992h . The standard Mathis, Rumpl, Nordsieck (MRN; iMathis et al.lll977f ) 
power-law was used for the grain size distributions. The dust temperature on the inner shell 
boundary and the optical depth were varied assuming an inverse square law for the spherical 
density distribution. The shell was assumed to extend to 10,000 times its inner radius. In 
order to model the dust shells of RSGs, we surveyed the numerical space of 3 parameters: i.e., 
the effective temperature of the central heating source (T e g), dust temperature on the inner 
shell boundary (T d ), and optical depth at 2.2 /im (t 22 ). The ranges and steps of parameter 
searches were 1500 K-4000 K and 100 K, 100 K-1800 K and 100K, and 0.001-0.05 and 0.001 
for T fj, Td and r 2 . 2 , respectively. Additionally, we inspected a couple of models with large 
optical depths (r 2 .2 > 0.05), but we confirmed, in such a range of large optical depths, that 
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there are no chances to fit the observational SEDs due to the strong absorption of silicate 
at 9.7 yum. We adopted the fits for which the sum of squares of the deviations between the 
observed and modeled fluxes (after scaling) were minimum. The mass-loss rate calculated 
by DUSTY posses general scaling properties in terms of the luminosity L, gas-to-dust mass 
ratio r g d and dust grain bulk density p s , assuming L = 10 4 L , r g d = 200 and p s = 3 g cm" 3 . 
The mass-loss rate is scaled in proport ion to L 3/,4 (r g dp s ) 1 ^ 2 (This quantity has roughly 30% 



inherent uncertainty; llvezic et al.lll999l ). To obtain the mass-loss rate of each RSG in Wdl, 
firstly we calculated the absolute luminosity of each RSG: we integrated flux densities by 
approximating the area under SED profiles by trapezoids, and the absolute luminosity was 
derived by assuming the distance of 3.5 kpc. The derived absolute luminosities (logL) are 
5.37, 6.03, 5.10 and 4.83 for W 237, W 26, W 20 and W 75, respectively. Then, using these 
absolute luminosities, we scaled the mass- loss rate, assuming r gd = 200 and p s = 3 g cm -3 . 
Consequently, the mass-loss rates (logM) are calculated to be —5.00, —4.65, —5.07 and 
—5.69 for W 237, W 26, W 20 and W 75, respectively. The best-fit input parameters of the 
DUSTY modeling, obtained mass-loss rates and absolute luminosities used for scaling are 
summarized in Table 7. 

Even though the derived mass-loss rates and lum inosities are within a typical range 



of the mass-loss rate of RSGs (IMauron fc Josselinll201ll ). W 75, which is a non-detection in 



both the SiO and H2O maser lines, exhibits the smallest luminosity and mass-loss rate, while 
two detections (W 237 and W 26) clearly exhibit larger absolute luminosities and mass-loss 
rates than W 75. According to the results of above analysis, one may expect that RSGs 
with maser emission are found predominantly in RSGs with relatively large mass-loss rates 
and absolute luminosities. 



5.2. Consideration using Previous Maser Observations 

As discussed in Section 5.1, the detection rates of maser emission of RSGs seem to be 
related to infrared properties. However, the number of RSGs in Wdl is unfortunately not 
sufficient enough to conclude something. Therefore, in this section, we extend our analysis 
to RSGs associated with other clusters (i.e., RSGC1, RSGC2, RSGC2 SW, Per OBI, Mc8). 
Since the galactocentric distance of these 5 clusters are, more or less, similar (except for 
Per OBI; see Table 6), we can expect that the effects of metallicity gradient in the galactic 



disk (IHenry fc Wortheylll999l ) are minimized. Unfortunately, at this moment, the distance 
information is based on radial velocities (i.e., distances are estimated by assuming the flat 
rotation model of the Galaxy). Therefore, the distances could include relatively large un- 
certainty [relative error of more than 50% according to the recent results of trigonometric 
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parallax measurements; Imai, H., in private communication; theoretically, however, the dis- 
tances to massive clusters may be improved in the future, because there are independent 



meth odologies to measure the distances to clusters (see, e.g.. |Perryman et al.lll998l ; lAn et al. 
20071 )]. Here, as a preliminary analysis, anyway we discuss based on the kinematic distances 
estimated from the radial velocities. 

In the five massive clusters mentioned above, in total 53 RSGs have been ide ntified until 
now ([Humphreys! 1 1 9 70l ; iFiger et al.ll2006l ; iDavies et al.ll2007l ; iDeguchi et al.ll2010l ). As well as 
we did in Section 5.1, we collected infrared photometric data of RSGs from DENIS, 2MASS 
PSC, Spitzer/ 'GLIMPSE, MSX, WISE All Sky Source Catalog and AKARI/IRC PSC. The 
collected infrared photometric data and SED plots are summarized again in Appendix B. 
Similarly to the analysis on Wd l, we adopted extinction corrections based on t he extinction 



coefficients given in literature (llndebetouw et al.l l2005t iLevesque et al.l 120051 ; IDavies et al. 



20071 . l2008l ). As we mentioned in Section 5.1, some authors recently use a slightly different 
formula for the conversion of the extinction values: we used Ak = 1.8 x Eh-k, while some 



authors recently use Ak = 1.5 x Eh-k (see, e.g., iNishiyama et al.l 120061 ) . We numerically 
evaluated the difference between the results obtained by the two formulas: for the majority 
of the samples with A K < 2.5 (this A K is based on the values given in Table Bl), the relative 
difference in the finally derived luminosities is less than 20%. For exceptional sources with a 
large extinction (Ak > 3.0), the relative difference occasionally goes beyond 40%, but such 
objects are quite rare (only two). Therefore, the selection of the conversion formula does not 
affect the results of our discussion. In Figure 4, we present the distribution of 29 RSGs in 
6 clusters (including Wdl) on the two-color diagram of log(F 2 i/F 12 ) and [K-12.13] colors. 
Although 5 RSGs in Per OBI (the red open diamonds) exhibit relatively small [K-12.13] 
values, there seems to be no other significant inhom ogeneity in the di stribution [since RSGs 
in Per OBI have been identified mainly in optical ( IHumphreysl ll97Cf ) . the sample could be 
somewhat biased]. Then, absolute luminosities and mass- loss rates were estimated using the 
same methodology as we adopted for Wdl. The obtained absolute luminosity and the best- 
fit parameters of DUSTY models are summarized in Table 8. Then, the results of previous 
maser observations of 53 RSGs in 5 clusters are summarized in Tables 9 and 10. 

Figure 5 shows the two-color diagrams of RSGs in 5 clusters and Wdl. The data points 
in Figure 5 are classified into either detections (red filled circles) or non-detections (blue open 
circles) in each maser line. In these two-color diagrams, we clearly see that the distributions 
of maser detections and non-detections are different in colors: i.e., detections are distributed 
in the upper-right of each panel, while the non-detections are distributed in lower-left of each 
panel. Since generally log(F 21 /F 12 ) and [K-12.13] are related to dust temperature and mid - 
infrared optical depth of the envelopes, respectively (see, e.g.. lLe Sidaner fc Le Bertrelll996l ). 
the concentration of maser detections suggests that maser emission is selectively detected 
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toward RSGs with a low dust temperature and large optical depth. On the other hand, we 
cannot see clear difference between the SiO, H 2 and OH maser lines in the distributions 
of detected and non-detected sources (even though the number of OH maser observations is 
not sufficient enough for comparison). 

Interestingly, however, one may point out a weak but possible difference between the 
SiO, H 2 and OH maser lines in the distributions of detections and non-detections in the 
logM versus logL plot given in Figure 6 (Note: in this figure we do not intend to discuss 
the correlation between the mass-loss rate and absolute luminosity; it is natural that there 
is a correlation between M and L in Figure 6, because both values are derived from the 
same SED. The main focus here is only the distributions of detections and non-detections in 
mass-loss rates and absolute luminosities). The detections in the SiO line distribute in a wide 
range of L, whereas the majority of the H 2 maser detections tend to locate in a region with 
large L (above, roughly 10 5 L ) except for one data point (St2-08) around logL ~ 4.55. Of 
course, we need to bear in mind that the distance includes a relatively large uncertainty, and 
the number of the samples is not large enough. However, as a possibility, one may interpret 
this difference for astrophysical reasons; for example, difference in excitation conditions of 
masers, such as, difference in excitation temperatures and critical densities. And if this 
difference is originated in the intrinsic astrophysical reasons, one may expect that maser 
lines are used to select a particular group of RSGs lying in a small range of mass-loss rate 
and luminosity. On the OH maser observation, at the moment we cannot see any tendency, 
because the number of observations are clearly too small to say something. Since the physical 
conditions required for pumping OH maser are also different from those for SiO and H 2 
masers, one may expect that there will be a recognizable difference from SiO and H 2 masers 
if the number of OH observations is increased. 

Using the photometric data collected, we inspected the distribution of detections and 
non-detections of ma se r searches on the Q1 -Q2 diagram, which has been suggested by 



Messineo et al.l (120121 ). iMessineo et al.l (120121 ) proposed a new source selection method for 



identifying mass-losing evolved stars, which is based on photometric data from 2MASS (J, 
H, K-band magnitudes) and GLIMPSE (8.0/mi-band magnitude) and on the Ql and Q2 
parameters. These two parameters are independent of interstellar extinction. The Ql pa- 
rameter [Ql = ( J — H) — 1.8 x (H — K)] is a measure of the deviation from the reddening 
vector in the J — H versus H — K plane. Early-type stars have Ql values around 1, K-giants 
around 0.4 mag, and dusty circumstellar envelopes of evolved stars genera te even smaller 



Ql v alues. On average, higher mass-loss rates produce smaller Ql values ( iMessineo et al. 



2012J). Similarly, the Q2 parameter [Q2 = (J -K) -2.69 x (K - [8.0]); here, [8.0] means the 
GLIMPSE 8 /zm-band magnitude] measures the deviation from the reddening vector in the 
plane J — K versus K — [8.0]. On average, higher mass-loss rates imply smaller Q2 values 



-14- 



as well as Ql ( IMessineo et all 120121 ) 



In Figure 7, we present the distribution of the present RS G samples on the Q1 -Q2 
diagram (red an d blue marks) ; the sam ples of Mira-like AGB stars (IMessineo et al.ll2002l ) and 
Galactic RSGs (jVerheyen et al.ll2012l ) are also plotted together (gray marks) for comparison. 
Unfortunately, only 13 data points of the present samples can be plotted on the diagram, 
because for the majority of stars the GLIMPSE 8 /xm data are not available due to source 
confusion (the present RSG samples are take n from star clus t ers, in which RSGs roll up in a 
small region). Nevertheless, as suggested by lVerheyen et all (120121 ). we see a clear tendency 
in the diagram: the detections in the SiO J = 1-0 lines exhibit negative Ql values, while 
non-detections exhibit relatively larger Ql value. Although we see the detection in the H2O 
maser line (St2-08) at (Ql, Q2) = (0.33,0.61), this presumably suggests that the SiO maser 
line is more sensiti ve to the Q-paramet ers than the H 2 maser lin e (St2-08 is nega t ive in the 
SiO maser search, iDeguchi et all 120101 ) . In the previous work by IVerheyen et al.l ( 120121 ) on 
Galactic RSG samples, the detections in the SiO J = 2-1 line are lying below Q2 = —1.9. 
However, in the present case we find one detection (Mc8-05) above Q2 = —1.9. Anyway, a 
larger statistics is required to confirm the suggested boundary on the Q1-Q2 diagram. 



6. Summary 

In this paper, we have reported the result of an ATCA observations of the Wdl region in 
the SiO v — 1, J — 1-0 and H 2 6i6-523 lines. We analyzed infrared SEDs of 57 RSGs in 6 
clusters including Wdl to investigate the correlations between maser and infrared properties. 
The main results of this research are summarized below: 

1. The SiO v = 1, J = 1-0 and H2O 6i6-523 lines are detected toward 2 of 4 known 
RSGs in Wdl, and the large velocity ranges of detected maser lines are consistent with 
the RSG status. The velocity ranges of the SiO maser line are wider than that of the 
H2O line; this fact may suggest that the RSG envelopes are deviated from a spherical 
symmetry. 

2. RSGs with maser emission seem to exhibit relatively large log(i 7 2 1 /F 12 ) and [K-12.13] 
colors compared to non-detections. The mass-loss rates derived from dust radiative 
transfer modeling suggests that RSGs with maser emission exhibit a relatively large 
mass-loss rate compared to RSGs with no maser emission. 

3. RSGs with SiO maser emission homogeneously distribute in L, whereas those with 
H2O maser emission tend to distribute in a region with large L values (above, roughly 
4 x 10 4 L ). 
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4. The distribution of the detections and non-detection on the Q I-Q2 diagram is roughly 
consistent with the previous study by IVerheyen et al.l ( 120121 ). even though one data 
point of a SiO maser detection is found above Q2 = —1.9. 
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A. New SiO Maser Detection toward IRAS 16105-4205 



We detected the SiO v = 1, J = 1-0 line toward IRAS 16105-4205 for the first time. 
The peak velocity and peak flux are —83 km s -1 and 8.8 Jy, respectively. The velocity- 
integrated flux of the d etection is 28. 2 Jy k m s" 1 . The SiO v = 1, J = 2-1 line has been 
detected previously by IHaikala et al.l (ll994[ ). In Figure Al, we present the line profile of 
the v = 1, J = 1-0 line superimposed on that of the v — 1, J — 2-1 line (IHaikala et al. 



1994T ). IRAS 16105 -4205 was previous ly classified as a carbon star based on near-infrared 
photometric colors (IFouque et al.lll992l ). However, since the object is detected in SiO maser 
lines, we suspect th at it may be an oxygen -rich star rather than a carbon star. In fact, 
Dollerv et all f)l987h and lchen fc Gaol J2002h respectively detected the OH 1612 MHz maser 
emission and 9.7 (im silicate feature toward this object, although a H 2 maser search was 
negative (jDeguchi et al.lll989[ ). 



B. Photometric Data of RSGs Associated with Massive Clusters 

Photometric data and extinction coefficients used for the analysis are summarized in 
Tables B1-B5. Spectral energy distribution diagrams and the results of DUSTY modeling 
are given in Figures B1-B3. 



- 16 - 



REFERENCES 

Alexander, M. J., Kobulnicky, H. A., Clemens, D. P., Jameson, K., Pinnick, A., & Pavel, M. 
2009, ApJ, 137, 4824 

An, D., Terndrup, D. M., Pinsonneault, M. H., Paulson, D. B., Hanson, R. B., & Stauffer, 
J. R. 2007, ApJ, 655, 233 

Asaki, Y., Deguchi, S., Imai, H., Hachisuka, K., Miyoshi, M., & Honma, M. 2010, ApJ, 721, 
267 

Boboltz, D. A., & Claussen, M. J. 2004, ApJ, 608, 480 
Borgman, J., Koornneef, J., & Slingerland, J. 1970, A&A, 4, 248 

Brandner, W., Clark, J. S., Stolte, A., Waters, R., Negueruela, I., & Goodwin, S. P. 2008, 
A&A, 478, 137 

Chen, P.-S., & Gao, Y.-F. 2002, Chinese J. Astron. Astrophys., 2, 169 

Cho, S.-H., Kaifu, N., & Ukita, N. 1996, Astron. Astrophys. Suppl. Ser., 115, 117 

Choi, Y. K., Hirota, T., Honma, M., Kobayashi, H., Bushimata, T., et al. 2008, PASJ, 60, 
1007 

Clark, J. S., Negueruela, I., Crowther, P. A., & Goodwin, S. P. 2005, A&A, 434, 949 

Clark, J. S., Negueruela, I., Davies, B., Larionov, V. M., Ritchie, B. W., Figer, D. F., 
Messineo, M., Crowther, P. A., & Arkharov, A. A. 2009, A&A, 498, 109 

Crowther, P. A., Hadfield, L. J., Clark, J. S., Negueruela, I., Vacca, W. D., Krabbe, A., 
Genzel, R., Eckart, A., Najarro, F., Lutz, D., Cameron, M., Kroker, H., Tacconi- 
Garman, L. E., Thatte, N., Weitzel, L., Drapatz, S., Geballe, T., Sternberg, A., & 
Kudritzki, R. 2006a, MNRAS, 372, 1407 

Crowther, P. A., Hadfield, L. J., Clark, J. S., Negueruela, I., & Vacca, W. D. 2008, MNRAS, 
385, 544 

Davies, B., Figer, D. F., Kudritzki, R.-P. et al. 2007, ApJ, 671, 781 
Davies, B., Figer, D. F., Law, C. J. et al. 2008, ApJ, 676, 1016 
Deguchi, S., Nakada, Y., Forster, J. R. 1989, MNRAS, 239, 825 
Deguchi, S., Nakashima, J., Zhang, Y., et al. 2010, PASJ, 62, 391 



-17- 



Deul, E., Epchtein, N., Borsenberger, J. 1995, Information & on-line data in astronomy, 203, 
73 

Dollery, M. E., Gaylard, M. J., Cohen, R. J. 1987, IAUS, 122, 215 

Egan M. P., Price S. D., Moshir, M. M., Cohen, M., & Tedesco, E., 1999, The Mid-Course 
Space Experiment Point Source Catalog Version 1.2 Explanatory Guide (AFRLVS- 
TR-1999-1522) 

Egan M.P, Price S.D., Kraemer K.E., Mizuno D.R., Carey S.J., Wright CO., Engelke C.W., 
Cohen M., & Gugliotti G. M. 2003, VizieR Online Data Catalog, 5114, 

Figer, D. F., MacKenty, J. W., Robberto, M. et al. 2006, ApJ, 643, 1166 

Fouque, P., Le Bertre, T., Epchtein, N., Guglielmo, F., Kerschbaum, F. 1992, A&AS, 93, 
151 

Haikala, L. K., Nyman, L.-A., & Forsstrom, V. 1994, Ap&SS, 103, 107 

Humphreys, R. M. 1970, ApJ, 160, 1149 

Henry, R. B. C. & Worthey, G. 1999, PASP, 111, 919 

Indebetouw, R., Mathis, J. S., Babler, B. L., Meade, M. R., Watson, C, Whitney, B. A., 
Wolff, M. J., Wolfire, M. G., Cohen, M., Bania, T. M., Benjamin, R. A., Clemens, D. 
P., Dickey, J. M., Jackson, J. M., Kobulnicky, H. A., Marston, A. P., Mercer, E. P., 
Stauffer, J. R., Stolovy, S. R., & Churchwell, E. 2005, ApJ, 619, 931 

Ishihara, D., Onaka, T., Kataza, H., Salama, A., Alfageme, C, Cassatella, A., Cox, N., 
Garcfa-Lario, P., Stephenson, C, Cohen, M., Fujishiro, N., Fujiwara, H., Hasegawa, 
S., Ita, Y., Kim, W., Matsuhara, H., Murakami, H., Miiller, T. G., Nakagawa, T., 
Ohyama, Y., Oyabu, S., Pyo, J., Sakon, I., Shibai, H., Takita, S., Tanabe, T., Uemizu, 
K., Ueno, M., Usui, F., Wada, T., Watarai, H., Yamamura, I., & Yamauchi, C. 2010, 
A&A, 514, Al 

Ivezic, Z., Nenkova, M., & Elitzur, M. 1999, User Manual for DUSTY, University of Kentucky 
Internal Report 

Jewell, P. R., Batrla, W., Walmsley, C. M., & Wilson, T. L. 1984, A&A, 130, LI 

Jiang, B. W., Deguchi, S., Yamamura, I., Nakada, Y., Cho, S. H., & Yamagata, T. 1996, 
ApJS, 106, 463 

Jiang, B. W., Deguchi, S., & Ramesh, B. 1999, PASJ, 51, 95 



Josselin, E., & Plez, B. 2005, ESO Astrophysics Symposia, 2005, 405 
Kamohara, R., Deguchi, S., Miyoshi, M., & Shen, Z.-Q. 2005, PASJ, 57, 341 
Kothes, R., & Dougherty, S. M. 2007, A&A, 468, 993 
Le Sidaner, P., & Le Bertre, T. 1996, A&A, 314, 896 
Levesque, E. M. 2010, New A Rev., 54, 1 

Levesque, E. M., Massey, P., Olsen, K. A., Plez, B., Josselin, E., Maeder, A., & Meynet, G. 
2005, ApJ, 628, 973 

Luna, A., Mayya, Y. D., Carrasco, L., Rodriguez- Merino, L. H., & Bronfman, L. 2009, 
RevMexAA, 37, 32 

Mathis, J. S., Rumpl, W. & Nordsieck, K. H. 1977, ApJ, 17, 25 

Mauron, N., & Josselin, E. 2011, A&A, 526, A156 

Mengel, S., & Tacconi-Garman, L. E. 2009, Ap&SS, 324, 321 

Mengel, S., & Tacconi-Garman, L. E. 2007, A&A, 466, 151 

Messineo, M., Habing, H. J., Sjouwerman, L. O., Omont, A., & Menten, K. M. 2002, A&A, 
393, 115 

Messineo, M., Menten, K. M., Churchwell, E., & Habing, H. 2012, A&A, 537, A10 
Nakashima, J., & Deguchi, S. 2006, ApJ, 647, L139 
Nakashima, J., & Deguchi, S. 2007, ApJ, 669, 446 

Negueruela, I., Gonzalez- Fernandez, C, Marco, A., Clark, J. S., & Martinez-Nunez, S. 2010, 
A&A, 513, A74 

Negueruela, I., Clark, J. S., & Ritchie, B. W. 2010, A&A, 516, A78 

Negueruela, I., Gonzalez-Fernandez, C, Marco, A., & Clark, J. S. 2011, A&A, 528, A59 

Nishiyama, S., Nagata, T., Kusakabe, N., Matsunaga, N., Naoi, T., et al. 2006, ApJ, 638, 
839 

Ossenkopf, V., Henning, Th., & Mathis, J. S. 1992, A&A, 261, 567 

Palagi, F., Cesaroni, R., Comoretto, G., Felli, M., & Natale, V. 1993, A&AS, 101, 153 



- 19 - 



Pashchenko, M. L, Rudnitskij, G. M., Samodurov, V. A., & Tolmachev, A. M. 2006, Astron. 
Astrophys. Trans., 25, 399 

Perryman, M. A. C, Brown, A. G. A., Lebreton, Y., Gomez, A., Turon, C, Cayrel de 
Strobel, G., Mermilliod, J. C, Robichon, N., Kovalevsky, J., & Crifo, F. 1998, A&A, 
331, 81 

Reid, M. J., Menten, K. M., Zheng, X. W., Brunthaler, A., Moscadelli, L., Xu, Y., Zhang, 

B. , Sato, M., Honma, M., Hirota, T., Hachisuka, K., Choi, Y. K., Moellenbrock, G. 

A. , & Bartkiewicz, A. 2009, ApJ, 700, 137 

Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, Astronomical Data Analysis Software 
and Systems IV, ASP Conference Series, eds, Shaw, R. A., Payne, H. E., & Hayes, J. 
J. E., 77, 433 

Scholz, M., & Wood, P. R. 2000, A&A, 362, 1065 

Skrutskie, M. F., Cutri, R. M., Stiening, R., Weinberg, M. D., Schneider, S., Carpenter, J. 
M., Beichman, C, Capps, R., Chester, T., Elias, J., Huchra, J., Liebert, J., Lonsdale, 

C, Monet, D. G., Price, S., Seitzer, P., Jarrett, T., Kirkpatrick, J. D., Gizis, J. E., 
Howard, E., Evans, T., Fowler, J., Fullmer, L., Hurt, R., Light, R., Kopan, E. L., 
Marsh, K. A., McCallon, H. L., Tam, R., Van Dyk, S., & Wheelock, S. 2006, AJ, 131, 
1163 

Stothers, R. B., & Chin, C-W. 1999, ApJ, 522, 960 

Su, J. B., Shen, Z.-Q., Chen, X., Yi, Jiyune, Jiang, D. R., & Yun, Y. J. 2012, ApJ, 754, 47 
Szymczak, M., Woolak, P., Gerard, E., & Richards, A. M. S. 2010, A&A, 524, A99 
Takaba, H., Iwata, T., Miyaji, T., & Deguchi, S. 2001, PASJ, 53, 517 
Verheyen, L., Messineo, M., & Menten, K. M. 2012, A&A, 541, A36 

Werner, M. W., Roellig, T. L., Low, F. J., Rieke, G. H., Rieke, M., Hoffmann, W. F., Young, 
E., Houck, J. R., Brandl, B., Fazio, G. G., Hora, J. L., Gehrz, R. D., Helou, G., Soifer, 

B. T., Stauffer, J., Keene, J., Eisenhardt, P., Gallagher, D., Gautier, T. N., Irace, 
W., Lawrence, C. R., Simmons, L., Van Cleve, J. E., Jura, M., Wright, E. L., & 
Cruikshank, D. P. 2004, ApJS, 154, 1 

Wilson, W. J., & Barrett, A. H. 1972, A&A, 17, 385 



-20 - 



Wright, Edward L.; Eisenhardt, Peter R. M., Mainzer, A. K., Ressler, M. E., Cutri, R. M., 
Jarrett, T., Kirkpatrick, J. D., Padgett, D., McMillan, R. S., Skrutskie, M., Stanford, 
S. A., Cohen, M., Walker, R. G., Mather, J. C, Leisawitz, D., Gautier, T. N., Ill, 
McLean, I., Benford, D., Lonsdale, C. J., Blain, A., Mendez, B., Irace, W. R., Duval, 
V., Liu, F., Royer, D., Heinrichsen, I., Howard, J., Shannon, M., Kendall, M., Walsh, 
A. L., Larsen, M., Cardon, J. G., Schick, S., Schwalm, M., Abid, M., Fabinsky, B., 
Naes, L., & Tsai, C-W. 2010, AJ, 140, 1868 

Zhang, B., Reid, M. J., Menten, K. M., & Zheng, X. W. 2012a, ApJ, 744, 23 

Zhang, B., Reid, M. J., Menten, K. M., Zheng, X. W., & Brunthaler, A. 2012b, A&A, 544, 
A42 



This preprint was prepared with the A AS IAT^X macros v5.2. 



-21 - 




Fig. 1. — Positions of FOVs for the SiO and H2O maser observations at 43 GHz and 22 GHz, 
respectively, superimposed on the J, H and K-bands color composite 2MASS image (image 
size is 8.8' x 5.9'; north is up, east is left). The blue and green circles represent FOVs for the 
SiO and H 2 observations, respectively. The diameters of the blue and green circles (i.e., 
HPBW of an ATCA 22 m antenna) are 66" and 126", respectively. Each FOV is numbered 
according to Table 1. Known RSGs with no maser emission are labeled with object names in 
white color. RSGs detected in the SiO and H 2 maser lines are labeled with object names 
in yellow color. 
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Fig. 2. — Spectra of the SiO v = 1, J = 1-0 (upper panels) and H 2 6i6-523 (lower panels) 
maser lines. The left and right panels show the line profiles of W 26 and W 237, respectively. 
The arrows indicate weak components (> 5a) shifted from the primary peak (see, text). 
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Fig. 3. — Spectral energy distributions of 4 known RSGs in Well. The blue open circles 
are data points of original values without interstellar extinction corrections, and the red 
filled circles are data points with interstellar extinction corrections. The light-blue and 
light-red triangles represent the lower limits of fluxes, respectively, before and after adopting 
interstellar extinction correction, and the light-blue and light-red crosses represent the flux 
values with a low quality flag, respectively, before and after adopting interstellar extinction 
correction. The black curves are the results of model fitting with the DUSTY code (see, text 
for details). 
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Fig. 4. — Distribution of RSGs in 6 massive clusters on the infrared two-color diagram. 
The member RSGs of Wdl are labeled with the arrows (W 26 is not indicated due to 
a large uncertainty in its K-band flux). In the log(F 2 i/Fi 2 ) color, F 2 \ and F\ 2 represent 
the MSX 21.34 /j,m and 12.13 fim flu xes, respectively. The [K-12.13] color is defined by 
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.fT-band magnitude and the MSX 12.13 fim flux. The interstellar extinction corrections were 
adopted (see Sections 5.1 and 5.2). 
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Fig. 5. — Distributions of detections and non-detections of maser observations in the SiO 
(J = 1-0, v = 1 and 2), H 2 (22 GHz) and OH maser (1612 MHz) lines on the infrared 
two-color diagrams. The definitions of colors are the same as Figure 4. The red filled and 
blue open circles represent detections and non-detections. The member RSGs of Wdl are 
indicated by the arrows. W 26 is not indicated in this plot due to large uncertainty in its 
photometric value. 
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Fig. 7. — Distributions of detections and non-detections of maser observations on the 
extinction-free Ql a nd Q2 parameters overlaid on a selected region of Figure 3 in 
Verheyen et al.l (120121 ) . The red filled circles represent RSGs detected only in the SiO maser 
line (and both H 2 and OH maser searches are negative). The red filled triangle represents 
RSG detected only in the H 2 maser line (and both SiO and OH maser searches are neg- 
ative). The open blue circles represent the RSGs that are negative in all three maser lines 
(SiO, H 2 and OH). The gray filled and empty hexa gons represent RSGs with and without 
maser detections, respectively ( Verheven et al. 20121). T he small dots represent SiO maser 
detections of Mira-like AGB stars ( Messineo et al. boo2h . The vertical dotted line represent 
the Q2 value of —1.9 (see, text). RSGs detected in maser lines are indicated by arrows along 
with object names. 
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Fig. Al — Spectrum of the SiO v = 1, J = 1-0 line of IRAS 16105-420 5 (thick line) 
super imposed on the spectrum of the SiO v = 1, J = 2-1 line (gray line; lHaikala et al. 
1994J ) . The v = 1, J = 1-0 line is detected for the first time. 
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Fig. Bl. — Spectral energy distributions of 53 known RSGs in 5 massive clusters. The 
notations of data points and model curves are the same as Figure 3. The names of parent 
clusters (small font) and object names (large font) are given in the upper-right corners of each 
panel. The names of objects with maser detections are indicated in boldface, and detected 
maser molecules are presented in the parenthesis. 
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Continue of Figure B2. 
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Table 1. Observational parameters 



Position of 


Center position 


Integration time 


Beam 


size* 


PA.* 


observation 


R.A. (J2000.0) 


Dec. (J2000.0) 


(hr) 










SiO Observation 


Position 1 


16:47:07.65 


-45:49:39.3 


0.65 


5.2" 


X 


2.5" 


5.6° 


Position 2 


16:47:04.20 


-45:50:38.3 


0.99 


5.1" 


X 


2.6" 


2.7° 


Position 3 


16:47:03.62 


-45:51:26.3 


0.99 


5.0" 


X 


2.6" 


4.7° 


Position 4 


16:47:05.34 


-45:52:24.3 


0.50 


8.1" 


X 


2.2" 


12.5° 


Position 5 


16:46:56.82 


-45:53:33.2 


0.50 


5.9" 


X 


2.4" 


14.7° 


Position 6 


16:47:26.12 


-45:52:17.3 


0.50 


8.5" 


X 


2.2" 


5.3° 


Position 7 


16:47:24.31 


-45:53:29.3 


0.50 


8.3" 


X 


2.2" 


9.3° 


Position 8 


16:46:43.06 


-45:49:53.0 


0.50 


5.8" 


X 


2.4" 


16.6° 


H2O Observation 


Position 1 


16:47:03.90 


-45:51:51.4 


2.03 


0.9" 


X 


0.4" 


-40.6° 


Position 2 


16:47:07.20 


-45:50:17.5 


0.58 


1.7" 


X 


0.3" 


-33.7° 



'Parameters of synthesized beams. 
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Table 2. List of known RSGs in Westerlund 1 



Source 


2MASS 


2MASS position 


name 


name 


R.A. (J2000.0) 


Dec. (J2000.0) 


W 237 


J16470309-4552189 


16:47:03.09 


-45:52:18.9 


W 26 


J16470540-4550367 


16:47:05.54 


-45:50:36.9 


W 20 


J16470468-4551238 


16:47:04.69 


-45:51:23.9 


W 75 


J16470892-4549585 


16:47:08.93 


-45:49:58.6 
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Table 3. Parameters of detected maser lines 



Source 


RA* 
(J2000.0) 


Dec* 
(J2000.0) 


-^pcak 

(Jy) 


hut 

(Jy km s _1 ) 


^pcak 

(km s- 1 ) 


rms 

(Jy) 


SiO (y = 


1, J = 1-0) 












W 237 
W 26 


16:47:03.07 
16:47:05.15 


-45:52:17.8 
-45:50:35.5 


1.54 
0.12 


16.82 
3.05 


-41.2 
-22.9 


0.027 
0.018 


H 2 (6ie 


-5 23 ) 












W 237 
W 26 


16:47:03.12 
16:47:05.41 


-45:52:19.2 
-45:50:36.9 


0.72 
0.73 


9.86 
3.29 


-51.0 
-39.3 


0.028 
0.094 



* Radio position obtained by fitting two-dimensional Gaussian function. 
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Table 4. Noise levels toward non-detected RSGs 



Source 


rms (SiO) 


rms (H 2 0) 




(Jy) 


(Jy) 


W 20 


0.005 


0.028 


W 75 


0.009 


0.094 
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Table 5. Summary of velocity measurements of Westerlund 1 



Measurement Averaged velocity 


Reference 


(km s- 1 ) 





SiO maser 


-40 


see Section 4 


H2O maser 


-44 


see Section 4 


CO bandhead 


-51 


Mensel & Tacconi-Garman (20091 


CO emission 


-48 ± 14 


Luna et al. (2009) 


HI absorption 


dJ -26 


Kothes & Douehertv (2007) 
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Table 6. Parameters of clusters 



Cluster 


Adopted dis. 
(kpc) 


Galactocentric dis. 
(kpc) 


Averaged VJ sr 
(km/s) 


Ref. 


Kin. dis. 
kot07* 
(kpc) 


Kin. dis. 
rei09+ 
(kpc) 


Averaged 
(mag) 


Wdl 


3.5 


4.5 


-39.9 


i 


3.0 


3.0 


1.33 


RSGC1 


6.6 


3.3 


120.9 


2 


6.9 


6.0 


2.57 


RSGC2 


5.5 


3.6 


96.6 


3 


5.2 


5.1 


1.40 


RSGC2 SW 


5.5 


.3.(5 


96.6 


3 


5.2 


5.2 


2.52 


Per OBI 


2.4 


9.5 


-38.5 


4 


3.3 


3.0 


0.31 


Mc8 


5.2 


3.3 


95.0 


3 


5.2 


5.2 


1.11 



* Calculated based on the Galaxy model of iKothes fe Dougherty! ||2007|) 
t Calculated based on the Galaxy model of iReid et al.l ll2009l ) 



References. — 1 see Section 4; • foavies et al ] |200S| ): foeguchi et al ] |20ich : ^Asaki et ail fcoid ) 
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Table 7. DUSTY model parameters and absolute luminosities 



Source 




T d 


T'2.2 


log M 


logL 


name 


(K) 


(K) 




(M© yr- 1 ) 


(L ) 


W 237 


3600 


600 


0.004 


-5.00 


5.37 


W 26 


3700 


500 


0.002 


-4.65 


6.03 


W 20 


3500 


500 


0.005 


-5.07 


5.10 


W 75 


3600 


1100 


0.004 


-5.69 


4.83 
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Table 8. DUSTY model parameters and absolute luminosities 



Source 


T c ft 


T d 


T2.2 


log M 


logL 


name 


(K) 


(K) 




(M© yr- 1 ) 


(L ) 


RSGC1 


F01 


3300 


300 


0.005 


-4.65 


5.41 


F02 


3600 


300 


0.004 


-4.58 


5.56 


F03 


3200 


400 


0.002 


-5.24 


5.22 


F04 


3900 


500 


0.001 


-5.46 


5.24 


F05 


3600 


400 


0.001 


-5.42 


5.22 


F06 


3300 


400 


0.001 


-5.61 


5.01 


F07 


3600 


1800 


0.001 


-6.36 


4.89 


F08 


3700 


300 


0.001 


-5.22 


5.30 


F09 


3100 


500 


0.003 


-5.30 


5.10 


F10 


3500 


200 


0.001 


-5.20 


5.09 


Fll 


3600 


1200 


0.002 


-5.73 


5.14 


F12 


3200 


1000 


0.001 


-6.04 


4.98 


F13 


4100 


500 


0.001 


-5.32 


5.40 


F14 


3500 


1500 


0.001 


-6.44 


4.67 


RSGC2 


St2-01 


3800 


1400 


0.003 


-6.06 


4.59 


St2-02 


3500 


1500 


0.001 


-6.47 


4.62 


St2-03 


3700 


900 


0.003 


-5.39 


5.20 


St2-04 


3700 


800 


0.001 


-5.91 


4.93 


St2-05 


3400 


1200 


0.001 


-6.33 


4.67 


RSGC2 SW 


St2-08 


3300 


600 


0.001 


-6.12 


4.55 


St2-ll 


3700 


600 


0.002 


-5.41 


5.12 


St2-14 


3500 


800 


0.001 


-5.90 


4.96 


St2-18 


3200 


500 


0.003 


-4.87 


5.64 


St2-26 


3300 


500 


0.001 


-6.14 


4.42 


Per OBI 


W Per 


3400 


600 


0.004 


-5.24 


5.09 


S Per 


3600 


900 


0.010 


-4.95 


5.27 


T Per 


3500 


1100 


0.001 


-6.31 


4.62 


V605 Cas 


3600 


500 


0.001 


-6.00 


4.56 


V778 Cas 


3500 


1200 


0.001 


-6.36 


4.61 


PR Per 


3300 


500 


0.001 


-6.01 


4.59 


FZ Per 


3400 


1300 


0.003 


-6.15 


4.48 


BD +59 372 


3100 


1500 


0.001 


-6.93 


4.08 


XX Per 


3700 


900 


0.002 


-5.73 


4.93 


HD 236947 


3400 


1500 


0.001 


-6.69 


4.35 


KK Per 


3500 


1300 


0.002 


-6.00 


4.85 


V550 Per 


3500 


1500 


0.002 


-6.51 


4.27 


PP Per 


3700 


1500 


0.001 


-6.67 


4.32 


BU Per 


3300 


400 


0.002 


-5.51 


4.83 


AD Per 


3400 


900 


0.001 


-6.01 


4.90 


SU Per 


3700 


700 


0.001 


-5.72 


5.10 


RS Per 


3500 


900 


0.003 


-5.53 


5.04 
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Table 8 — Continued 



Source 


T c tt 


T d 


T2.2 


log M 


logL 


name 


(K) 


(K) 




(M© yr- 1 ) 


(L ) 


V439 Per 


3600 


1500 


0.001 


-6.52 


4.54 


V403 Per 


3500 


1500 


0.001 


-6.70 


4.31 


V441 Per 


3200 


1000 


0.002 


-5.95 


4.78 


YZ Per 


3400 


800 


0.002 


-5.64 


4.84 


GP Cas 


3300 


800 


0.002 


-5.84 


4.78 


V648 Cas 


3700 


700 


0.001 


-5.79 


5.01 


Mc8 


Mc8-01 


3200 


600 


0.004 


-5.66 


4.56 


Mc8-02 


3200 


600 


0.001 


-6.27 


4.37 


Mc8-03 


3300 


900 


0.010 


-5.85 


4.12 


Mc8-04 


3300 


1400 


0.015 


-5.77 


4.35 


Mc8-05 


3200 


900 


0.002 


-6.18 


4.41 


Mc8-06 


3100 


1500 


0.003 


-6.16 


4.62 
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Table 9. Summary of previous maser observations of other clusters 1 



Source IRAS Coordinate Maser observation 



name 




name 


RA 


Dec 


SiO 


H 2 Q 


OH 








(J 2000.0) 


(J 2000.0) 








RSGC1 


F01 






18:37:56.30 


-06:52:32.2 


y 1 


n 1 


n 2 


F02 






18:37:55.28 


-06:52:48.3 


y 1 


n 1 


n 2 


F03 






18:37:59.73 


-06:53:49.2 


n 1 


n 1 


n 2 


F04 






18:37:50.88 


-06:53:38.1 


y 1 


n 1 


n 2 


F05 






18:37:55.52 


-06:52:11.6 


n 1 


n 1 


n 2 


F06 






18:37:57.45 


-06:53:25.3 


n 1 


n 1 


n 2 


F07 






18:37:54.31 


-06:52:34.5 


n 1 


n 1 


n 2 


F08 






18:37:55.19 


-06:52:10.7 






n 2 


F09 






18:37:57.76 


-06:52:22.1 


n 1 


n 1 


n 2 


F10 






18:37:59.53 


-06:53:31.9 


n 1 


n 1 


n 2 


Fll 






18:37:51.72 


-06:51:49.7 


n 1 


n 1 


n 2 


F12 






18:38:03.30 


-06:52:45.1 


n 1 


n 1 


n 2 


F13 






18:37:58.91 


-06:52:32.1 


y 1 


y 1 


y 2 


F14 






18:37:47.65 


-06:53:02.1 


n 1 


n 1 


n 2 


RSGC2 


St2-01 






18:39:19.88 


-06:01:48.0 


n 3 


n 3 




St2-02 






18:39:18.25 


-06:02:14.2 


n 3 


n 3 




St2-03 


IRAS 18366-0603 


18:39:19.59 


-06:00:41.6 


y 3 


y 3 


n 4 


St2-04 






18:39:24.60 


-06:02:13.9 


n 3 


n 3 




St2-05 






18:39:14.70 


-06:01:36.5 


n 3 


n 3 




RSGC2 SW 


St2-08 






18:39:07.75 


-06:03:20.2 


n 3 


y 3 




St2-ll 






18:39:05.59 


-06:04:26.6 


y 3 


y 3 




St2-14 






18:39:08.04 


-06:05:24.2 


n 3 


n 3 




St2-18 






18:39:02.38 


-06:05:10.6 


y 3 


y 3 




St2-26 






18:38:56.98 


-06:06:45.7 


y 3 


n 3 




Per OBI 


W Per 


IRAS 


02469+5646 


02:50:37.91 


56:59:00.1 


y 5 


n« 




S Per 


IRAS 


02192+5821 


02:22:51.73 


58:35:11.2 


y 7 


y 8 


y 9 


T Per 


IRAS 


02157+5843 


02:19:21.87 


58:57:40.3 


n io 






V605 Cas 


IRAS 


02167+5926 


02:20:22.48 


59:40:16.8 


n io 






V778 Cas 


IRAS 


01550+5901 


01:58:28.90 


59:16:08.7 


rilO 






PR Per 


IRAS 


02181+5738 


02:21:42.41 


57:51:46.0 


n io 






FZ Per 


IRAS 


02174+5655 


02:20:59.65 


57:09:30.0 


n io 






BD +59 372 


IRAS 


01561+6000 


01:59:39.67 


60:15:01.9 








XX Per 


IRAS 


01597+5459 


02:03:09.36 


55:13:56.6 




n 6 




HD 236947 


IRAS 


02036+5832 


02:07:12.06 


58:47:15.9 








KK Per 


IRAS 


02068+5619 


02:10:15.79 


56:33:32.7 








V550 Per 


IRAS 


02116+5754 


02:15:13.30 


58:08:32.3 








PP Per 


IRAS 


02135+5817 


02:17:08.23 


58:31:47.0 








BU Per 


IRAS 


02153+5711 


02:18:53.30 


57:25:16.8 




n 6 




AD Per 


IRAS 


02169+5645 


02:20:29.00 


56:59:35.2 








SU Per 


IRAS 


02185+5622 


02:22:06.89 


56:36:14.9 









Table 9 — Continued 



Source IRAS Coordinate Maser observation 



name 


name 


RA 


Dec 


SiO 


H 2 OH 






(J 2000.0) 


(J 2000.0) 






RS Per 


IRAS 02188+5652 


02:22:24.30 


57:06:34.4 




n 6 


V439 Per 


IRAS 02196+5658 


02:23:11.03 


57:11:58.3 






V403 Per 




02:23:24.11 


57:12:43.1 






V441 Per 


IRAS 02217+5712 


02:25:21.86 


57:26:14.1 






YZ Per 


IRAS 02347+5649 


02:38:25.42 


57:02:46.2 




n 6 


GP Cas 


IRAS 02360+5922 


02:39:50.44 


59:35:51.3 






V648 Cas 


IRAS 02473+5738 


02:51:03.95 


57:51:19.9 






Mc8 


Mc8-01 


IRAS 18258-1058 


18:28:35.32 


-10:56:36.6 


y 3 


n 3 


Mc8-02 




18:28:39.13 


-10:55:52.7 


n 3 


n 3 


Mc8-03 




18:28:52.08 


-10:57:57.7 


y 3 


n 3 


Mc8-04 




18:28:54.04 


-10:56:45.5 


n 3 


n 3 


Mc8-05 




18:28:58.46 


-10:56:07.0 


y 3 


n 3 


Mc8-06 




18:29:03.05 


-10:53:15.4 


n 3 


n 3 



References, 
i 20 lcM Wilson k, Barrettf 
i2001):f~ 
1 19991) 



1 Nakashima & Dei 



Cho et alj il996l ); fPaTagi et 



sguchil IpOOeh^Davies et al.1 j200Sl ^ toeguchi et al 
il972Mjiang et al.1 ( 199i 

t alj Il9! 



. . , ; e Takaba et al 
19931) : %zvmczak et al.1 j2010l) ; ^ Jiang et al 
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Table 10. Summary of previous maser observations of other clusters 2 



Cluster 


Line 


Detection/No. of RSGs observed 


Telescope & Ref. 


RSGCl 


SiO v = 1, J = 1-0 


4/13 


NR045m 1 




SiO v = 2, J = 1-0 


1/13 


NR045m 1 




H a O (6 16 -5 23 ) 


1/14 


NR045m 1 


OH 


2 n 3/2 J = 3/2 F = 1-2 


1/14 


VLA 2 


RSGC2 


SiO v = 1, J = 1-0 


1/5 


NR045m 3 




SiO v = 2, J = 1-0 


0/5 


NR045m 3 




H 2 (6i 6 -5 23 ) 


1/5 


NR045m 3 


OH 


2 n 3/2 J = 3/2 F = 1-2 


0/1 


NRA043m 4 


RSGC2 SW 


SiO v = 1, J = 1-0 


3/5 


NR045m 3 




SiO v = 2, J = 1-0 


2/5 


NR045m 3 




H 2 (616-523) 


3/5 


NR045m 3 


Per OBI 


SiO v = 1, J = 1-0 


2/7 


NR045m 5 - 6 ' 7 




SiO v = 2, J = 1-0 


1/7 


NR045m 5 - 6 ' 7 




SiO v = 3, J = 1-0 


0/1 


NR045m 6 




H 2 (6i 6 -5 23 ) 


1/7 


Medicina32m 8 ,KSRC34m 9 


OH 


2 n 3/2 J = 3/2 F = 1-2 


1/1 


NRT 10 


OH 


2 n 3/2 J = 3/2 F = 2-2 


1/1 


NRT 10 


Mc8 


SiO v = 1, J = 1-0 


3/6 


NR045m 3 




SiO ti = 2, J = 1-0 


3/6 


NR045m 3 




H 2 (616-533 ) 


0/6 


NR045m 3 


NR045m - 45m radio telescope, Nobeyama Radio Observatory; 




VLA - The Very Lar 


ge Array, National Radio 


Astronomy Observatory; 





Medicina32m - Medicina 32m Antenna, Istituto di Radioastronomia; 
KSRC34m - Kashima 34m Antenna, Kashima Space Research Center; 
NRT - Nangay radio telescope, Observatoire de Paris; 

NRA043m - 43-m radiotelescope, National Radio Astronomy Observatory 



Ref er ences. — ^Nakashima fc Deguc hi feOOrjl^Davies et al.l j2008h;toeguchi et al.l JioToh A Wilson fc Barrett 
1972[) ; 4jiang et al.l <199d); e Cho et al.1 <1996fl : 1jiang et al.l i ll999l)^PaIagi et al.l jl993h ; ^Takaba et al. 
200j) ; 1 HS2vmczak et al.l feOlcK ; 
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Table Bl. Infrared Photometric Data 1 



Source 


A K 


DENIS" 






2MASS 6 






name 


2.17 (im 


I 




J 


0J 


H 




K 


OK 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


Wdl 


W 237 


0.96 


10.25 


0.03 


5.08: 


0.02 


3.01 


0.27 


2.18 


0.28 


W 26 


1.33 


9.74 


0.04 


4.31 


0.52 


1.35 


1.32 


6.21f 




W 20 


1.55 


13.05 


0.03 


6.38 


0.02 


4.38: 


0.08 


2.61 


0.41 


W 75 


1.47 


14.03 


0.03 


6.93 


0.02 


4.67: 


0.02 


3.28 


0.27 


RSGC1 


F01 


2.58 






9.75 


0.03 


6.59 


0.02 


4.96 


0.02 


F02 


2.83 






9.90 


0.03 


6.70 


0.02 


5.03 


0.02 


F03 


2.46 






9.95 


0.03 


6.92 


0.02 


5.33 


0.02 


F04 


2.46 






9.66 


0.03 


6.80 


0.02 


5.34 


0.02 


F05 


2.77 






10.55 


0.03 


7.18 


0.03 


5.54 


0.03 


F06 


2.19 






9.87 


0.03 


7.04 


0.03 


5.61 


0.03 


F07 


2.33 






9.94 


0.03 


7.07 


0.03 


5.63 


0.02 


F08 


2.84 






10.77 


0.03 


7.33 


0.02 


5.65 


0.02 


F09 


2.44 






10.26 


0.03 


7.24 


0.02 


5.67 


0.02 


F10 


2.45 


18.47 


0.29 


10.18 


0.03 


7.22 


0.03 


5.71 


0.02 


Fll 


2.63 






10.47 


0.03 


7.33 


0.03 


5.72 


0.02 


F12 


2.57 






10.14 


0.03 


7.24 


0.02 


5.86 


0.03 


F13 


3.19 






10.91 


0.03 


7.72 


0.02 


5.96 


0.02 


F14 


2.29 






10.50 


0.03 


7.58 


0.02 


6.17 


0.02 


RSGC2 


St2-01 


1.45 


13.92 


0.03 


7.82 


0.02 


6.02 


0.03 


5.11 


0.02 


St2-02 


1.42 


14.35 


0.03 


8.35 


0.03 


6.21 


0.03 


5.26 


0.02 


St2-03 


1.39 


12.66 


0.03 


6.90 


0.02 


5.05 


0.02 


4.12 


0.27 


St2-04 


1.34 


12.85 


0.03 


7.27 


0.02 


5.46 


0.03 


4.50: 


0.04 


St2-05 


1.42 


13.94 


0.03 


8.22 


0.02 


6.21 


0.02 


5.24 


0.02 


RSGC2 SW 


St2-08 


1.64 


15.19 


0.05 


8.57 


0.02 


6.31 


0.02 


5.23 


0.02 


St2-ll 


4.58 






14.23 


0.04 


9.92 


0.02 


7.32 


0.02 


St2-14 


1.99 


16.03 


0.07 


8.53 


0.02 


6.05 


0.03 


4.82 


0.02 


St2-18 


1.87 


14.56 


0.04 


7.15 


0.03 


4.70 


0.02 


2.90 


0.19 


St2-26 


2.52 






12.03 


0.04 


8.52 


0.03 


6.95f 




Per OBI 


W Per 


0.45 






3.10 


0.20 


2.00 


0.17 


1.57 


0.25 


S Per 


0.47 






2.95 


0.28 


1.85 


0.24 


1.12 


0.22 


T Per 


0.35 






4.03 


0.25 


3.02 


0.17 


2.58 


0.24 


V605 Cas 


0.35 






3.95 


0.26 


3.04 


0.20 


2.48 


0.27 


V778 Cas 


0.19 






4.03 


0.30 


2.97 


0.21 


2.49 


0.32 


PR Per 


0.28 






3.56 


0.31 


2.68 


0.19 


2.25 


0.30 


FZ Per 


0.22 






3.80 


0.25 


2.91 


0.20 


2.48 


0.24 


BD+59 372 


0.28 






5.33 


0.25 


4.20 


0.23 


3.74 


0.26 


XX Per 


0.31 






3.44 


0.26 


2.48 


0.23 


1.97 


0.24 


HD 236947 


0.31 






4.02 


0.23 


3.13 


0.19 


2.62 


0.24 


KK Per 


0.31 






3.00 


0.20 


2.14 


0.17 


1.68 


0.21 
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Table Bl — Continued 



Source 


A K 


DENIS" 






2MASS 6 






name 


2.17 






J 




H 




K 






(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


V550 Per 


0.31 






4.72 


0.27 


3.66 


0.20 


3.19 


0.32 


PP Per 


0.31 






4.57 


0.25 


3.52 


0.22 


2.95 


0.24 


BU Per 


0.36 






3.68 


0.29 


2.68 


0.20 


2.19 


0.23 


AD Per 


0.31 






3.38 


0.24 


2.48 


0.21 


1.94 


0.26 


SU Per 


0.23 






2.82 


0.25 


1.93 


0.18 


1.46 


0.22 


RS Per 


0.29 






3.05 


0.22 


2.11 


0.19 


1.56 


0.21 


V439 Per 


0.21 






4.13 


0.29 


3.22 


0.23 


2.69 


0.30 


V403 Per 


0.24 






4.48 


0.27 


3.46 


0.22 


2.99 


0.29 


V441 Per 


0.28 






3.47 


0.26 


2.47 


0.22 


2.04 


0.24 


YZ Per 


0.26 






3.26 


0.26 


2.30 


0.23 


1.91 


0.26 


GP Cas 


0.31 






3.69 


0.27 


2.44 


0.21 


1.94 


0.24 


V648 Cas 


0.45 






3.83 


0.34 


2.71 


0.22 


2.11 


0.27 


Mc8 


Mc8-01 


1.05 


16.37 


0.09 


7.90 


0.03 


5.98 


0.03 


4.88 


0.02 


Mc8-02 


0.88 


13.41 


0.05 


7.74 


0.02 


5.86 


0.04 


4.97 


0.02 


Mc8-03 


1.40 






9.78 


0.02 


7.54 


0.04 


6.23 


0.02 


Mc8-04 


1.51 






9.30 


0.02 


6.90 


0.04 


5.59 


0.02 


Mc8-05 


1.24 


17.44 


0.14 


8.88 


0.04 


6.80 


0.04 


5.57 


0.02 


Mc8-06 


0.59 


11.06 


0.02 


6.84 


0.03 


5.26 


0.03 


4.49 


0.02 



a For saturating sources, the measured fluxes are calibrated by PSF-fitting. 

b For saturating sources, the measured fluxes are calibrated by 1-D radial profiles. 

: The colon represents uncertain detection. 

tThe dagger represents 95% confidence upper limit. 
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Table B2. Infrared Photometric Data 2 



Source WISE" 

name 3.4 /^m cr 3 . 4 ^ m 4.6 (im (J4.6/im 12 /^m CTi2^m 22 /im o"22^tm 

(mag) (mag) (mag) (mag) (mag) (mag) (mag) (mag) 



RSGC1 



F01 


3.37 


0.25 


2.49 


0.16 


-0.71 


0.16 


-2.10 


0.02 


F02 


3.88 


0.09 






-1.22 


0.07 


-2.77 


0.01 


F03 


4.93 


0.13 


4.20 


0.11 


1.50 


0.07 


-0.50 


0.07 


F04 


5.54 


0.08 


4.72 


0.05 


1.27 


0.06 


-0.83 


0.02 


F05 


3.89 


0.12 


4.05 


0.02 


0.99 


0.07 


-0.36 


0.02 


F06 


5.41 


0.06 


4.42 


0.04 


2.24 


0.02 


-0.17 


0.02 


F07 


5.57 


0.11 


5.26 


0.45 


3.78 


0.42 


4.89f 




F08 


3.89 


0.12 


4.05 


0.02 


0.99 


0.07 


-0.36 


0.02 


F09 


















F10 


4.68 


0.09 


3.86 


0.06 


2.14 


0.03 


-0.04 


0.02 


Fll 


5.49 


0.06 


4.42 


0.05 


2.50 


0.03 


0.82 


0.04 


F12 


4.92 


0.06 


3.97 


0.04 


2.55 


0.02 


1.31 


0.02 


F13 


3.47 


0.16 


2.02 


0.08 






-1.98 


0.01 


F14 


5.15 


0.05 


4.63 


0.04 


3.43 


0.02 


2.26 


0.03 


RSGC2 


St2-01 


4.96 


0.07 


4.22 


0.05 


2.98 


0.02 


1.56 


0.03 


St2-02 


4.81 


0.07 


4.22 


0.04 


3.14 


0.03 


1.61 


0.03 


St2-03 


2.85 


0.18 


1.93 


0.11 


-0.26 


0.03 


-1.26 


0.02 


St2-04 


4.51 


0.09 


3.91 


0.05 


1.81 


0.02 


0.55 


0.03 


St2-05 


4.62 


0.08 


4.11 


0.05 


2.85 


0.02 


1.61 


0.03 


RSGC2 SW 


St2-08 


4.82 


0.06 


4.19 


0.05 


3.12 


0.02 


1.45 


0.09 


St2-ll 


6.33 


0.05 


4.07 


0.06 


1.33 


0.07 


-0.63 


0.03 


St2-14 


5.00 


0.05 


4.37 


0.04 


2.27 


0.02 


0.66 


0.04 


St2-18 


2.28 


0.25 


1.44 


0.01 


-1.27 


0.06 


-2.12 


0.02 


St2-26 


5.72 


0.05 


4.69 


0.04 


3.25 


0.02 


1.76 


0.04 


Per OBI 


W Per 


2.56 


0.03 


1.70 


0.01 


-1.73 


0.21 


-2.89 


0.001 


S Per 










-3.11 


0.29 


-4.32 


0.001 


T Per 


2.43 


0.12 


1.85 


0.07 


1.01 


0.01 


-0.19 


0.01 


V605 Cas 


2.21 


0.32 


1.58 


0.09 


0.71 


0.02 


-0.44 


0.02 


V778 Cas 


2.33 


0.06 


1.94 


0.01 


0.70 


0.01 


0.32 


0.01 


PR Per 


2.01 


0.29 


1.52 


0.04 


0.49 


0.01 


-0.56 


0.02 


FZ Per 


2.27 


0.16 


1.49 


0.12 


0.62 


0.01 


0.07 


0.01 


BD+59 372 


3.19 


0.13 


2.62 


0.05 


2.69 


0.01 


2.60 


0.02 


XX Per 






1.41 


0.05 


-0.94 


0.04 


-1.86 


0.01 


HD 236947 


2.66 


0.12 


1.86 


0.09 


1.99 


0.01 


1.68 


0.02 


KK Per 


2.12 


0.27 


1.68 


0.07 


-0.27 


0.01 


-0.79 


0.01 


V550 Per 


2.94 


0.14 


2.41 


0.09 


1.90 


0.01 


1.33 


0.01 


PP Per 


2.60 


0.26 


2.08 


0.12 


2.01 


0.02 


1.64 


0.03 


BU Per 


1.94 


0.48 


1.30 


0.20 


-0.64 


0.02 


-1.99 


0.004 


AD Per 


1.88 


0.47 


1.30 


0.19 


-0.17 


0.02 


-0.91 


0.01 


SU Per 






1.38 


0.21 


-1.29 


0.02 


-1.88 


0.01 



Table B2— Continued 



Source WISE" 



name 


3.4 (im 0-3. 4 ^ m 
(mag) (mag) 


4.6 fim 

(mag) 


C4.6^tm 
(mag) 


12 fim 
(mag) 


012/jm 
(mag) 


22 /im 

(mag) 


o"22/im 
(mag) 


RS Per 










-1.23 


0.03 


-2.43 


0.01 


V439 Per 


2.56 


0.26 


2.03 


0.07 


1.13 


0.01 


1.01 


0.01 


V403 Per 










1.94 


0.02 


1.78 


0.02 


V441 Per 






1.33 


0.05 


0.06 


0.01 


-0.95 


0.01 


YZ Per 


2.26 


0.01 


1.81 


0.02 


-1.03 


0.01 


-1.72 


0.01 


GP Cas 










-0.34 


0.02 


-1.30 


0.01 


V648 Cas 


2.02 


0.34 






-0.70 


0.01 


-1.79 


0.004 


Mc8 


Mc8-01 


3.95 


0.08 


2.83 


0.06 


1.52 


0.02 


0.48 


0.02 


Mc8-02 


5.14 


0.06 


4.39 


0.04 


3.31 


0.02 


2.28 


0.03 


Mc8-03 


5.53 


0.06 


4.64 


0.03 


2.82 


0.02 


1.31 


0.02 


Mc8-04 


4.84 


0.07 


3.24 


0.05 


1.50 


0.02 


0.49 


0.03 


Mc8-05 


5.31 


0.06 


4.16 


0.04 


2.57 


0.02 


1.31 


0.03 


Mc8-06 


3.62 


0.09 


3.29 


0.05 


2.35 


0.02 


0.73 


0.02 



a For saturating sources, the measured fluxes are calibrated by non-saturated wings of their profiles. 
^The dagger represents 95% confidence upper limit. 
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Table B3. Infrared Photometric Data 3 



Source GLIMPSE" 



name 


3.6 /Ltm 


"3.6fim 


4.5 /urn 


"4.5fim 


5.8 /im 


05.8/im 


8.0 Atm 


"S.Ofim 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


RSGC1 


F01 


















F02 






3.54 


0.07 


2.69 


0.05 






F03 


4.23 


0.11 






3.64 


0.04 






F04 


4.21 


0.04 


4.13 


0.05 


3.60 


0.04 






F05 


5.24 


0.22 






3.85 


0.03 


5.05 


0.24 


F06 










3.83 


0.03 






F07 






4.62 


0.07 


3.86 


0.03 


3.25 


0.04 


F08 










3.89 


0.03 






F09 


4.47 


0.05 


4.55 


0.07 


3.80 


0.03 






F10 


4.70 


0.03 


4.64 


0.04 


3.96 


0.03 






Fll 


4.43 


0.10 






3.89 


0.03 






F12 










4.17 


0.03 


3.99 


0.04 


F13 


















F14 


5.89 


0.12 


5.88 


0.12 


4.65 


0.03 


4.36 


0.03 


RSGC2 


St2-01 


4.65 


0.08 


4.74 


0.07 


4.34 


0.03 


4.26 


0.03 


St2-02 


4.56 


0.05 






4.23 


0.03 


4.11 


0.03 


St2-03 


















St2-04 


















St2-05 


4.72 


0.19 






4.23 


0.03 


4.02 


0.03 


RSGC2 SW 


St2-08 


4.56 


0.06 






4.22 


0.02 


4.22 


0.03 


St2-ll 


















St2-14 


3.89 


0.08 


3.98 


0.05 


3.61 


0.03 






St2-18 


















St2-26 


















Mc8 


Mc8-01 


4.02 


0.21 


4.95 


0.24 


3.30 


0.04 






Mc8-02 


4.65 


0.06 






4.27 


0.03 


4.21 


0.02 


Mc8-03 






5.91 


0.17 


4.46 


0.03 


3.98 


0.02 


Mc8-04 


4.43 


0.14 






3.58 


0.03 


3.27 


0.04 


Mc8-05 


5.03 


0.26 


4.85 


0.07 


4.35 


0.03 


4.22 


0.03 


Mc8-06 


3.99 


0.06 


4.31 


0.05 


3.83 


0.03 


3.54 


0.03 



a For saturating sources, the measured fluxes are calibrated by non-saturated wings of their profiles. The 
aperture corrections of point sources were made using the effective correction factors given by Reach et 
al. (2005). 



Table B4. Infrared Photometric Data 4 



Source MSX 

name 4.29 /im 0-4.29/jm 4.35 /im 0-4.35/im 8.28 ;itm o-g.28^tm 12.13 fim (Ji2.l3/im 14.65 /^m cri4.65^m 21.34 /im (T2i.34/im 

(Jy) (Jy) (Jy) (Jy) (Jy) (Jy) (Jy) (Jy) (Jy) (Jy) (Jy) (Jy) 



Wdl 



W 237 


43.40 


4.21 


45.50 


4.19 


61.17 


2.51 


100.60 


5.03 


82.95 


5.06 


101.39 


6.08 


W 26 


143.27 


12.18 


171.67 


15.11 


216.37 


8.87 


326.91 


16.35 


254.14 


15.50 


240.82 


14.45 


W 20 


30.33 


3.25 


37.82 


3.56 


37.06 


1.52 


57.67 


2.88 


52.18 


3.18 


62.56 


3.75 


W 75 






19.05 


2.04 


13.27 


0.54 


14.43 


0.72 


11.25 


0.69 


11.35 


0.68 


RSGC1 


FOl 
F02 
F03 






6.19 


0.59 


5.77 


0.24 


11.91 


0.60 


9.68 


0.59 


9.82 


0.59 


F04 
F05 
F06 
F07 






5.96 


0.57 


5.42 


0.22 


10.71 


0.54 


8.75 


0.53 


6.99 


0.42 






2.98 


0.34 


3.21 


0.13 


4.98 


0.25 


3.45 


0.21 


3.19 


0.19 


F08 
F09 
F10 
Fll 






5.90 


0.57 


5.14 


0.21 


7.82 


0.39 


6.70 


0.41 


8.74 


0.52 






3.49 


0.38 


2.66 


0.11 


3.82 


0.19 


2.79 


0.17 


1.77 


0.11 


F12 






3.47 


0.38 


1.87 


0.08 


2.49 


0.12 


1.44 


0.09 


1.15 


0.07 


F13 


16.28: 


7.65 


8.69 


0.81 


8.71 


0.36 


14.17 


0.71 


12.30 


0.75 


14.17 


0.85 


F14 






2.28: 


0.30 


0.97 


0.04 


0.84 


0.05 


0.64 


0.04 


0.41: 


0.03 


RSGC2 


St2-01 










2.03 


0.08 


1.93 


0.11 


1.21 


0.08 






St2-02 










1.63 


0.07 


1.48 


0.09 


0.67: 


0.06 






St2-03 






15.74: 


2.27 


13.14 


0.54 


19.91 


1.00 


15.38 


0.94 


13.19 


0.79 


St2-04 


18.16: 


9.03 


10.07: 


1.54 


5.12 


0.21 


6.27 


0.31 


4.02 


0.25 


2.49 


0.17 


St2-05 










1.71 


0.07 


1.52 


0.10 


1.08 


0.08 






RSGC2 SW 


St2-08 










1.38 


0.06 


1.54 


0.10 


1.01 


0.07 






St2-ll 










7.35 


0.30 


16.61 


0.83 


16.56 


1.01 


14.25 


0.86 


St2-14 










3.38 


0.14 


4.88 


0.25 


3.92 


0.24 


2.94 


0.19 


St2-18 


27.77 


2.83 


48.62 


5.01 


31.24 


1.28 


56.86 


2.84 


45.89 


2.80 


35.74 


2.14 


St2-26 










1.37 


0.06 


1.83 


0.11 


1.34 


0.09 







Per OBI 



Table B4— Continued 



Source 














MSX 












name 


4.29 /urn 

(Jy) 


""4.29/tm 

(Jy) 


4.35 fim 

(Jy) 


0-4.35^111 
(Jy) 


8.28 /im 

(Jy) 


<78.28/im 

(Jy) 


12.13 fan 

(Jy) 


CT i2.l3/im 

(Jy) 


14.65 fim 

(Jy) 


C"14.65/im 

(Jy) 


21.34 fim 

(Jy) 


CT 21.34/im 

(Jy) 





W Per 


50.99 


5.15 


63.11 


5.87 


63.45 


2.60 


79.47 


3.97 


62.61 


3.82 


68.76 


4.13 


S Per 


147.27 


12.81 


186.16 


16.57 


162.82 


6.68 


274.56 


13.73 


185.62 


11.32 


178.67 


10.72 


T Per 


29.18: 


3.73 


19.93 


2.13 


10.57 


0.43 


8.66 


0.43 


6.57 


0.40 


4.75 


0.30 


V605 Cas 


16.99: 


2.46 


27.32 


2.71 


12.54 


0.51 


11.32 


0.57 


7.69 


0.47 


7.93 


0.48 


V778 Cas 


- 


- 


32.23 


3.35 


10.22 


0.42 


7.20 


0.38 


4.94 


0.32 


2.89: 


0.51 


PR Per 


26.75: 


3.56 


28.78: 


3.28 


14.36 


0.59 


11.78 


0.59 


8.90 


0.54 


8.26 


0.51 


FZ Per 


- 


- 


- 


- 


12.19 


0.50 


10.90 


0.54 


6.52 


0.40 


4.13 


0.27 


BD+59 372 










2.80 


0.11 


1.57 


0.10 


0.89: 


0.07 






XX Per 


























HD 236947 


























KK Per 


























V550 Per 






16.67: 


2.32 


5.19 


0.21 


3.51 


0.22 


1.81 


0.15 






PP Per 


32.92: 


9.81 


12.56: 


2.10 


5.48 


0.22 


2.74 


0.17 


2.31 


0.16 






BU Per 


33.14 


3.38 


34.91 


3.56 


33.10 


1.36 


36.70 


1.83 


26.30 


1.60 


30.21 


1.81 


AD Per 






35.80 


3.37 


20.75 


0.85 


20.44 


1.02 


13.12 


0.80 


11.25 


0.69 


SU Per 






20.58 


2.41 


43.77 


1.79 


40.03 


2.00 


24.11 


1.47 


30.08 


1.80 


RS Per 






73.08 


6.72 


57.59 


2.36 


59.47 


2.97 


42.78 


2.61 


41.09 


2.47 


V439 Per 






17.79: 


2.38 


7.33 


0.30 


5.18 


0.27 


3.20 


0.20 


2.39 


0.17 


V403 Per 


30.91 


3.21 


15.85: 


2.25 


5.31 


0.22 


2.71 


0.16 


1.83 


0.12 






V441 Per 


39.42 


3.86 


44.46 


4.31 


19.32 


0.79 


16.81 


0.84 


13.26 


0.81 


11.10 


0.68 


YZ Per 


37.91 


4.25 


42.70 


4.18 


33.62 


1.38 


31.78 


1.59 


21.23 


1.29 


20.18 


1.21 


GP Cas 


45.68 


4.29 


48.57 


4.57 


22.97 


0.94 


23.51 


1.18 


17.14 


1.05 


16.10 


0.97 


V648 Cas 


54.72 


5.03 


50.40 


4.74 


31.87 


1.31 


34.81 


1.74 


26.89 


1.64 


24.57 


1.47 



Mc8 



Mc8-01 


3.80: 


0.48 






4.66 


0.19 


5.61 


0.28 


4.68 


0.29 


3.67 


0.22 


Mc8-02 


2.82: 


0.42 






1.28 


0.05 


1.04 


0.05 


0.79 


0.05 


0.73 


0.05 


Mc8-03 


10.78: 


2.12 






3.39 


0.14 


4.12 


0.21 


3.12 


0.19 


3.24 


0.19 


Mc8-04 


3.23: 


0.45 






3.21 


0.13 


3.64 


0.18 


2.80 


0.17 


2.04 


0.12 


Mc8-05 










1.44 


0.06 


2.17 


0.11 


1.63 


0.10 


1.02 


0.07 


Mc8-06 


5.00 


0.55 






2.46 


0.10 


2.27 


0.11 


1.45 


0.09 


0.75 


0.05 



: The colon represents uncertain detection. 
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Table B5. Infrared Photometric Data 5 



Source AKARI 



name 


9.0 /urn 


cg^tm 


18.0 /im 


o"i8^tm 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


Wdl 


W 237 


91.65 


7.35 


113.80 


0. 


67 


W 26 


296.50 


7.68 


449.70 


4. 


31 


W 20 


61.41 


0.28 


111.70 


0. 


.67 


W 75 


24.26 


1.64 








RSGC1 


F01 


11.47 


0.34 








F02 


12.69 


0.17 








F03 


6.46 


0.54 








F04 


5.75 


0.03 








F05 


7.86 


0.45 








F06 


5.17 


0.58 








F07 












F08 


7.86 


0.45 








F09 












F10 












Fll 


4.01 


0.70 








F12 


2.16 


0.12 








F13 


13.30 


0.41 








F14 


0.91 


0.01 








RSGC2 


St2-01 


2.64 


0.13 


1.52 


0. 


08 


St2-02 


1.51 


0.05 


0.56 


0. 


05 


St2-03 


14.16 


0.26 


14.24 





.17 


St2-04 


6.98 


1.43 


3.68 


0. 


08 


St2-05 


1.45 


0.01 


0.66 


0. 


.09 


RSGC2 SW 


St2-08 


1.30 


0.16 








St2-ll 


8.62 


0.41 


12.31 


0. 


24 


St2-14 






2.80 





.14 


St2-18 


33.98 


0.41 


36.78 





.08 


St2-26 


1.45 


0.20 


1.17 


0. 


09 


Per OBI 


W Per 


77.36 


0.82 


73.29 


0. 


77 


S Per 


349.70 


3.68 


285.50 


11. 


80 


T Per 


9.80 


0.12 


5.77 


0. 


09 


V605 Cas 


12.00 


0.14 


7.64 


0. 


17 


V778 Cas 


9.61 


0.01 


3.85 


0. 


04 


PR Per 


13.01 


0.08 


8.09 





10 


FZ Per 


12.03 


0.12 


4.98 





.05 


BD+59 372 


2.52 


0.02 


0.66 





.02 


XX Per 


64.90 


1.57 


28.47 


0. 


40 


HD 236947 


4.69 


0.02 


1.35 


0. 


02 


KK Per 


19.18 


0.08 


9.50 


0. 


.08 
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Table B5 — Continued 



Source AKARI 

name 9.0 [tm ag/xm 18.0 [tm <Ti8//m 

(Jy) (Jy) (Jy) (Jy) 



V550 Per 


4.63 


0.03 


1.67 


0.01 


PP Per 


4.76 


0.05 


1.43 


0.01 


BU Per 


36.59 


0.51 


30.27 


0.51 


AD Per 


20.81 


0.15 


11.90 


0.13 


SU Per 


45.48 


0.47 


30.69 


0.41 


RS Per 


66.55 


0.88 


47.17 


0.85 


V439 Per 


6.93 


0.04 


2.29 


0.02 


V403 Per 


4.81 


0.02 


1.33 


0.03 


V441 Per 






11.91 


0.15 


YZ Per 


37.24 


0.36 


25.05 


0.08 


GP Cas 


24.89 


0.40 


15.95 


0.04 


V648 Cas 


31.81 


1.29 


24.37 


0.30 


Mc8 


Mc8-01 


5.17 


0.11 






Mc8-02 


1.43 


0.05 


0.67 


0.05 


Mc8-03 






2.62 


0.13 


Mc8-04 






3.00 


0.58 


Mc8-05 






1.17 


0.04 


Mc8-06 











